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Abstract 



A theory of spin polarized tunneling spectroscopy based on a scattering theory 
is given for tunneling junctions between ferromagnets and d-wave supercon- 
ductors. The spin filtering effect of an exchange field in the insulator is also 
treated. We clarify that the properties of the Andreev reflection are largely 
modified due to a presence of an exchange field in the ferromagnets, and 
consequently the Andreev reflected quasiparticle shows an evanescent-wave 
behavior depending on the injection angle of the quasiparticle. Conductance 
formulas for the spin current as well as the charge current are given as a func- 
tion of the applied voltage and the spin-polarization in the ferromagnet for 
arbitrary barrier heights. It is shown that the surface bound states do not 
contribute to the spin current and that the zero-bias conductance peak ex- 
pected for a d-wave superconductor splits into two peaks under the influence 
of the exchange interaction in the insulator. 



PACS Numbers : 74.50.+r, 74.72.-h, 74.80.Fp 
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I. INTRODUCTION 



The transport properties in hybrid structures between ferromagnets and superconduc- 
tors have received considerable theoretical and experimental attentions. Interest in such 
structures includes spin-dependent spectroscopy of superconductors and possible device ap- 
plications. Since the Cooper pairs in spin singlet superconductors are formed between up 
and down spins, the high density of spin injection through a tunneling barrier induces a 
spin imbalance. This non-equilibrium state is expected to result in a suppression of the 
critical temperature and the critical current density in the superconductor. A large num- 
ber of experimental studies on spin-polarized tunneling have already been performed using 
conventional metal superconductors such as Al and Nb about 20 years agou. However, the 
recent discovery of so-called colossal magneto-resistance (CMR) in Mn oxides compound has 
aroused new interest in this fieldii, because hybrid structure fabrication of the spin-polarized 
ferromagnets with high-T^ superconductors is now possible using these materials^'!. 

On the other hand, the properties of ferromagnet/insulator/superconductor (F/I/S) 
and ferromagnet /ferromagnetic- insulator/superconductor (F/FI/S) junctions have been an- 
alyzed based on the assumption that the conductance spectra correspond to the density 
of states (DOS) of the superconductor weighted by the spin polarizationBi'i. A theory for 
F/I/S junctions based on a scattering method has been presented by de Jong and Beenekeri, 
and new aspects of Andreev reflection have been revealecL and also detailed comparisons 
between theory and experiments have been accomplishedBll2l. However, these results are 
restricted to isotropic s-wave superconductors. 

In contrast to s-wave superconductor cases, at the interface of a (ii.2_y2-wave supercon- 
ductor, zero-energy states fZES) are formed due to the interference effect of the internal 
phase of the pair potentialtll. Tunneling theory for (ia,2_j,2-wave superconductors has al- 
ready been presented by extending the BTK formulali^ to include the anisotropy of the pair 
potentiatJcJ. The theory predicts the existence of zero-bias conductance peak (ZBCP) 
which reflects the formation of the surface bound states on the d-wave superconductors. In 
this paper, an exchange interaction is introduced on the normal side of the junction and 
on the insulator in order to analyze the spin polarized tunneling effects. The bound-state 
condition and tunneling spectroscopy of ferromagnet /(i- wave superconductor junctions have 
already been analyzed in two papers0'0. They have revealed several important features in 
charge transport. Here we will argue that the properties of the Andreev reflectionEl is largely 
modifled due to the presence of the exchange interaction. In particular, the existence of an 
evanescent type of the Andreev reflection, which is referred to as virtual Andreev reflection 
(VAR), is explained for the flrst time (see Ref. 27). This process has signiflcant roles on the 
transports especially for junctions between half-metallic ferromagnets and superconductors. 
The conductance formulas for the charge and the spin currents are presented based on the 
scattering method by fully taking account of the VAR process. The merit of a formula 
based on the scattering methods is that the conductance spectra can easily be calculated 
for arbitrary barrier heights cases without the restriction of the high-barrier limit. The spin 
current is, we believe, the most important physical quantity in spin injection devices based 
on the following two reasons: one is that the spin current gives a direct criteria to estimate 
the effect of the spin imbalance induced by the tunneling current, the other is that the 
charge and the spin conductivity may illuminate the study of electron systems that undergo 
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spin-charge separation, such as Tomonaga-Luttinger hquids and possibly underdoped high- 
Tc superconductorsli3iHiJ. We will also analyze ferromagnetic insulator effects, which includes 
the spin-filtering effectill, due to the presence of an exchange field in the insulator. It is 
shown that a spin- dependent energy shift during the tunneling process induces a splitting of 
the ZBCP. Based on the detailed analysis of the conductance spectra, we propose a simple 
method to distinguish the broken time-reversal symmetry (BTRS) states inducement at the 
surfac(ilB from spin-dependent tunneling effects. The implications of the ferromagnetic 
insulator effects on tunneling experiments of high-T^ superconductors and a proposal for 
possible device applications are also presented. 



II. FORMULATION 



For the model of formulation, a planar F/FJ/S" junction with semi-infinite electrodes in 
the clean limit is assumed. A flat interface is assumed to be located at x = 0, and the insu- 
lator for up [down] spin is described by a potential V|^](a;) m[|](a;) = (Vq — [+]Ub)S{x)}, 
where S{x), Vq and Ub are the 5-function, a genuine barrier amplitude and an exchange 
amplitude in the barrier, respectively. The effective mass m in the ferromagnet and in the 
superconductor are assumed to be equal. For the model of the ferromagnet, we adopt the 
Stoner model where the effect of the spin polarization is described bvthe one-electron Hamil- 
tonian with an exchange interaction similarly to the case of Ref.ErllSllZl. For the description 
of the dx2-y2-waye superconductor, we apply the quasi-classical approximation where the 
Fermi energy Ep in the superconductor is much larger than the pair potential following the 
model by Brudei§§. The effective Hamiltonian (Bogoliubov-de Gennes equation) is give 
by 



Hoix) - pU{x) 
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Here, E is the energy of the quasiparticle, U{x) is the exchange potential given by UQ{—x) 
{U > 0) where Q{x) is the Heaviside step function, p is 1 [-1] for up [down] spins, A{x,9) 
is the pair potential and Hq{x) = —h'^V'^/2m + V{x) — Ep. To describe the Fermi surface 
difference in F and S, we assume Ep = EpN for a; < and Ep = Eps for a; > 0. The pair 
potential A{x,9) is taken as A{9)Q{x) for simplicity. The number of up [down] spin electrons 
is described by The polarization and the wave- vector of quasiparticles in the 



ferromagnet for up [down] spin are expressed as P| 
and kj' 



Efn+U 



[Pi 



Efn-U] 



f;v,T =1 fciv,T \= ^wi^FN + U) [kM,i =1 k^,i \= ^^^{EpN-U)i respectively!. 
We assume the quasiparticle injection of up spin electrons at an angle 6]s to the interface 
normal as shown in Fig. |l|. Four possible trajectories exist; they are Andreev reflection 
(AR), normal reflection (NR), transmission to superconductor as electron-like quasiparticles 
(ELQ), and transmission as hole-like quasiparticles (HLQ). The spin direction is conserved 
for NR but not for AR. When the superconductor has (ij:2_y2-wave symmetry, the effective 
pair potentials for ELQ and HLQ are given by A+ = Aq cos 2(6*5—/?) and A_ = Aq cos 2(6'5'-|- 
/9), respectively, where (3 is the angle between a-axis of the crystal and the interface normal. 
Results for various pairing symmetries are obtained by setting proper values to A+ and A_ 
similarly to the previous formula^Ull. The wave vectors of ELQ and HLQ are approximated 
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by ks =1 ks |~ \f^^^ following the model by Andreevlii. Since translational symmetry 
holds along the y-axis direction, the momentum components of all trajectories are conserved 
(/cAf,t sin On = kN,i sin 9a = ks sin 9s)^. Note that 9n is not equal to 9a except when U = 0, 
which means retro-reflectivity of AR is broken. Such novel behavior is a consequence of the 
fact that in the presence of an exchange field the BCS paring is formed not strictly between 
states of equal but opposite fc-vectors, the so-called Fulde-Ferrell effectH. The wave-function 
in the ferromagnet {x < 0) for up [down] spin with injection angle 9n is described by 
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where the signs of the x-components of fc7v,t[i] aiid k' are the reversed of each other. 
The reflection probabilities of the two processes are obtained by solving Eq. (|ip and by 
connecting the wave-function and its derivative at a; = 0. 

Next, we will simply explain the Fermi surface effect by assuming up spin injection. 
Various kinds of reflection process are expected depending on the values of EpN, Es and 
U. For example, when fc^ < k^,^, total reflection (| b-i[i]{E,9N) P= 1) occurs when 9n > 
sin~^(/cs//c7v,t) = ^ci'^'^- In this case, the net currents of the spin and the charge from the 
ferromagnet to the superconductor vanish. On the other hand, when /cjv | < ks < kN,i, the 



x-component of wave-vector in AR process {yk%^^ — k^ sin 9s) becomes purely imaginary 

for 9ci > 9m > sin~^(/cAr^j^/fc7v,t) = ^c2- In this case, although transmitted quasiparticles 
from ferromagnet to superconductor do propagate, the Andreev reflected quasiparticles do 
not propagate (VAR process). A finite amplitude of the evanescent AR process still exists 
(I ai{E,9N) P> 0) and the net currents of the spin and the charge from the ferromagnet 
to the superconductor do not vanish. It is easy to check the conservation laws for the 
charge^the excitation, and the spin on the VAR process following the method presented 
in Ref.li3. The existence of the VAR process has not been treated in the one- dimensional 
modeli because it is a peculiar feature of a two or three dimensional F/S interface. 

The conductance of the junctions are obtained by extending previous formula to include 
the effect of spini'lE. In the following, consider a situation where kN,i < ks < k^,^- To 
analyze the transport properties of an F/J/S" junction, two kinds of conductance spectrum 
are introduced. The conductance for the charge current is defined by the charge flow induced 
by the up [down] spin quasiparticle injection and is given by 
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The conductance for the spin current is defined by the spin imbalance induced by the up 
[down] spin quasiparticle injection, 
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(for 0c2 <| 0w |< Oci) 
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(for 0a < I 07V I < vr/2) 
=0. 

The Andreev reflected quasiparticles positively contribute to the charge current, but since 
their spins are reversed, they have negative contribution to the spin current. Second terms 
in r.h.s. of Eqs. (|^) and (^ do not have finite contribution on net current in the VAR 
process, since the corresponding A2 is purely imaginary. The normalized total conductance 
spectra for the charge current OqiE) and the spin current OgiE) are given by 
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In the above, Rjy, ag^i^[ii{E) and crs^|[j](-E') correspond to the conductance when the super- 
conductor is in the normal state and the spin-resolved normalized conductance spectra for 
charge and spin, respectively. The net polarization Jp{eV) as a function of the bias voltage 
V is give by 

J. ... I-oodEa.{E){f{E-eV)-fm 
MeV) j^^dEa,{E){fiE - eV) - /(E)}' 

where f{E) is the Fermi distribution function. Since the convolution with f{E) gives only 
a smearing effect in the conductance spectra, the temperature is set to zero in the following 
discussions. 

In the above formulation, we have neglected the self-consistency of the pair potential 
in order to get analytical formulas^. However, the present formula is easily extended to 
include this effect simply by replacing T± with T±{x) \x=o, where T±{x) follows the Ricatti 
equations described by 

f+(x) = \-A+{x)Tlix) - Alix) + 2Ef+(x)l , (15) 
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Here the spatial dependence of the pair potential is assumed as A±(x) (functions of x). 

The most important differences in the present formula from previous ones are; i) a novel 
formula for the non-linear spin current, ii) a capability to treat the ferromagnetic insulator 
effects based on the scattering method, iii) the introduction of the breakdown in the retro- 
refiectivity of the AR process and consequently the vanishing of the propagating AR (VAR 
process) . In particular, the concept of the VAR process is a new physical process presented 
in this paper. If we would not accept the existence of this process, the total reflection 
independent of E is naively expected. Since finite transmission is possible in this angle 
region {6c2 <| 6* |< Od) above Tc, this total reflection would induce a sudden decrease of the 
conductance just below for highly polarized ferromagnets junctions. As far as we know, 
no trends for such effect has been reported thus far. This fact may be the direct evidence for 
the existence of the VAR process. The VAR process is shown to have an important role on 
the Josephson current in superconductor /ferromagnet /superconductor junctions, because 
the evanescent wave carries a net Josephson current in this configuratiorS. Note that the 
suppression mechanism of the AR process presented here is essentially different from that 
discussed in one-dimensional model where the contribution of the AR to the net current is 
simply governed by the ratio 
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III. RESULTS 



A. Effects of polarization 

In this subsection, to reveal the influence of the polarization on the tunneling conduc- 
tance spectra, we assume F/I/S junction by setting f/^ = (Zo,| = Zq^i = Zq). At first, let 
us discuss several analytical results obtained from above formulation in order to check the 
validity of the formula. When U = 0, the ferromagnet reduces to a normal metal, and as 
expected (Jq{E) reproduces the results of Ref. and (Js{E) vanishes. For half- metallic fer- 
romagnets {U = Epf^f), the Fermi-surface for the down spins has shrunk to zero. In this case, 
the VAR process occurs for all 9^- Under the condition of VAR, ag^E, 9n) = crs{E, 6^) ap- 
plies, which corresponds to the fact that the tunneling current is completely spin-polarized. 
Furthermore, the conductance spectra in the energy gap {E <| A_|_ |, <| A_ |) become 
completely zero [aq{E,0]^) = as{E,9iy) = 0]. In the tunneling limit {H — > oo) and in the 
absence of VAR, ag^i[i]{E , 9) gives the angle resolved surface DOS of an isolated supercon- 
ductor. Then (Tq{E) converges to the surface DOS weighted by the tunneling probability 
distributional. At this limit, we can reproduce a well-known result that the ratio of the_peak 
heights in the spin-resolved spectra directly reflect the polarization in the ferromagnetu. On 
the other hand, as^^[\]{E, 9) reduces to a function similar to the surface DOS, but where the 
divergence at the energy levels of the surface bound states is missing. 

Next, the calculated results based on above formula are presented for c?a.2_y2-wave super- 
conductors. In the following, we assume EpN = Eps- Figures]^ and ^show the conductance 
spectra of charge current for the transparent limit {Zq = 0, /? = 0) and high-barrier case 
(Zq = 5, /3 = 7r/4) as the function of exchange interaction X{= U/Efn)- For X = 0, 
results in Refill are reproduced. However, as X increases, the conductance inside the gap 
(I |< Aq) is largely reduced for both cases. Especially, the ZBCP disappears for the 
half-metallic ferromagnet case. Since the spin-polarization has such a drastic influence on 
the ZBCP, the height of ZBCP can be used in principle as a measurement of the magnitude 
of the spin polarization. Figure ^ shows the difference of the spin current and the charge 
current when X=0.85, = 5 and /3 = 7r/4. It is clear that the ZBCP is not present for the 
spin current. This corresponds to the fact that the charge current components corresponding 
to the ZES are carried by condensed Cooper pairs in the superconductor, and therefore they 
do not contribute to the spin imbalance. As a result, the spin current becomes relatively 
insensitive to the orientation of the junctions. Figure |^ shows the conductance spectra for 
the spin current as the function of spin polarization (Zq = 5). It is clear that the spin 
current increases as X becomes larger. Note that (Js{E) is larger than unity around E = Aq 
when X ^ 1. This corresponds to the fact that the peak in the DOS has an influence even 
for the spin current. 

Next, the net polarization Jp{eV) is calculated for d^2_y2-wsive superconductors as a 
function of the orientation (jS) when T = 0. Four lines of Fig. |^ show the results for various 
values of barrier parameter when eV = 2Aq. It is clear that the orientational effect is much 
smaller compared to the effect of Zq. In the same figure, results for s-wave superconductors 
(A+ = A_ = Aq independent of 9^) are also shown as closed dots. The large deviations 
of dx2-y2-'wave from s-wave for small values of Zq are originated from the distribution of 
the pair amplitude in A;-space. As the barrier parameter becomes larger, the spin injection 
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efficiency becomes to be insensitive to the symmetry of the pair potentiaL 



B. Spin filtering effects and the ZBCP splitting 

It has been experimentally verified that a ferromagnetic semiconductor used as the in- 
sulator in tunneling junctions works as a ferromagnetic barrier. Since the transmission 
probabilities for up and down spins are not equal, an spin-filtering effect is expected to be 
realizedi'0. Also it has been theoretically verified that a ferromagnetic insulator placed in the 
vicinity of superconductor induces a spin-splitting on the DOS of s-wave superconductors^. 
In the following, we will analyze the infiuence of the exchange interaction existing inside the 
insulator on the transport properties based on the formulation described in Sec. ||. 

Figure shows the response of the conductance spectra <Jg{E) on the exchange interaction 
in the insulator when X = 0. ZBCP splittings are obtained for finite exchange amplitude 
{Ub) cases. As Ub is increased and consequently as the difference between Zq -^ and Zq^^ 
becomes larger, the amplitude of the splitting becomes larger and the two peaks become 
broader and smaller. The peaks in the gap disappear when the difference between and 
Zo^l becomes prominent. To see more clearly these trends, the spin-resolved conductance 
spectra ag^i[i]{E) and <Js{E) for Zo_| = 2.5 and Zq^^ = 7.5 are plotted in Fig. |^. The spectra 
for up [down] spins are shifted for lower [higher] energy level. Furthermore, (Ts{E) becomes 
finite even though X = in the ferro magnet. In order to check the effect of the polarization. 
Fig. 1^ shows the response of the charge current as a function of polarization X for a fixed 
barrier parameter. The spin polarization in the ferromagnet induces the imbalance of the 
peak heights, thus the ratio of the splitted peak heights can be used as a criteria for the 
spin-polarization. 

These results are interpreted as follows: i) The peaks corresponding to the up [down] spin 
components are shifted because of the energvgain (loss) during the tunneling process, ii) 
Since this energy gain (loss) has /c-dependenceGil, the peak becomes broader comparing to the 
magnetic-field induced peak splitting (see below), iii) The amplitude of the peak splitting 
depends on the genuine barrier amplitude Vq as well as the exchange amplitude Ub- For 
example, the splitted peaks merge into a single peak at the tunneling limit (Vq — > cx)) even 
if Ub is kept constant, iv) The current corresponding to the ZBCP is carried by the Cooper 
pair in the superconductor as described in the previous subsection. This corresponds to the 
fact that the AR process is the second-lowest order tunneling process which requires both up 
and down spins tunneling. Hence, as Zq^^ becomes larger and as the tunneling probability for 
down spins are suppressed, the conductance peaks and the AR process are rapidly reduced 
even if Zo,| is kept zero, v) The spin current is increased as Ub is raised from zero even 
if X in the ferromagnet is kept at zero (unpolarized). This feature directly corresponds to 
the spin-filtering effect that the spin-selective tunneling occurs due to the presence of the 
exchange field in the insulator. 

Next, various types of the ZBCP splitting expected for d-wave superconductors and 
their polarization effects are analyzed. Mainly two possibilities other than the ferromag- 
netic insulator effects have been proposed for the origins of the ZBCP splitting on high- Tc 
superconductor junctions. One is the Zeeman effect due to an applied magnetic field, and 
the other is the inducement of the BTRS states such as dx2_y2+is-wa:Ve. The conductance 
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spectra in an applied magnetic field is calculated from above formula by simply using the 
relation 



(^qls]{E) = (^qis],^iE - flBH) + CTgisuiE + fisH), (17) 

where ^bH is the Zeeman energy. Calculated charge conductance spectra for (ia.2_y2-wave 
superconductor as a function of X are shown in Fig. [TO. The amplitude of the splitting is 
linear to the applied field independent of the barrier heights. Moreover, since the energy 
shift induced by the magnetic field does not have fc-dependence, the broadening of the peaks 
are not observed. The ratio of the splitted peak heights simply reflects the polarization 
in the ferromagnet, which is consistent with the results by Tedrow and Meservey0. On 
the other hand, (Jg{E) for (ij,.2_j^2+is-wave superconductor is calculated by setting A± = 
Aq cos 2(6^5 T P) + i^s- Calculated charge conductance spectra for various X values are 
shown in Fig. |rT[ The amplitude of the splitting is almost equivalent to the amplitude of 
the s-wave component. The shape of the spectrum without the polarization (X = 0) is quite 
similar to that shown in Fig. |I^ (X = 0). As X becomes larger, the heights of the two peaks 
are reduced, which is consistent with that shown in Fig. ^. On the other hand, differently 
from Figs. |^ and |TD|, since the peak splitting is not induced by spin-dependent effects in this 
case, the polarization in the ferromagnets does not yield an imbalance in the peak heights. 
Thus the heights of the two peaks are reduced symmetrically. 

The responses of the ZBCP on the variation of the polarization and the applied magnetic 
field are summarized as follows: i) the peak splitting due to the ferromagnetic insulator and 

the Zeeman effect are spin dependent. Therefore, the polarization in the ferromagnet 
induces the asymmetrical splitting of the ZBCP. ii) The amplitude of the peak splitting is 
linear to the applied field in the case of the Zeeman effect. However, it is non-linear in the 
cases of the ferromagnetic insulator effect^ and the BTRS statesi. In particular, the peak 
splittings are expected even in the absence of the applied field for these two cases, iii) The 
combination of the BTRS states and the Zeeman effect induces an additional peak splitting, 
that is, the ZBCP splits into four peaks. However, the combination of the Zeeman and the 
ferromagnetic insulator effects yields two peaks. 

The experimental observations of the ZBCP splitting have been reported for normal 
metal / high- Tc superconductor junctionsBBIlii! It is a really interesting experiment 
to observe the same features by using ferromagnet s/high-T^ superconductor junctions in 
order to distinguish the spin- dependent effects from the BTRS states inducement. Recently, 
Sawa, et.al. have detected an asymmetric magnetic field response in Lao.eySro.saMnOs / 
YBa2Cu307_5 junctions^. The qualitative features on the magnetic field responses of their 
junctions are consistent with F/FI/S with dx2-y2-waye explained here. Detailed comparison 
between above formulas and their experiments is strongly expected. 

Finally, a simple proposal is given for a possible device application utilizing the ferromag- 
netic insulator effects. The thickness of the insulator is the order of Inm in usual tunneling 
junctions. Since the controlling of properties in such a thin layer requires high technology, as 
far as we know, not so many experimental trials have been accomplished thus far. However, 
as shown in this paper, a small change in the insulator property causes a drastic change on 
the transport properties. Therefore, the controlling of the barrier properties is one of the 
most promising methods to create new functional devices. For example, consider a F/FI/S 
junction with a dx^-y^-wave superconductor {(3 = 7r/4). The sharp ZBCP is drastically 
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modified as the difference between and becomes larger as shown in Fig. |^. This means 
that, for a fixed bias voltage, a large response in current is expected due to a small variation 
in the exchange interaction in the insulator. This response is applicable for the high-sensitive 
magnetization measurement of a thin insulating film by inserting the film into a junction 
as a tunneling barrier. If the exchange interaction is sensitive to the external field, this 
effect can be used as a magnetic sensor. Alternatively, if the magnetization of the insulator 
shows a hysteresis on the external field variation, a memory function can be realized. The 
current gain of the junction as a function of the external field is largely enhanced by using 
a superconductor /ferromagnetic insulator /superconductor junction with (i-wave, because 
negative- conductance regions are expected just beside the ZBCP in this configurationili§. 
Differently from conventional superconducting memories based on a fiux-quantum logic, a 
large-scale integration circuit may be possible based on the present principle. 



IV. SUMMARY 



In this paper, the conductance spectra for the charge and the spin currents under the 
influence of the exchange interaction have been calculated based on the scattering method. 
The influence of the spin polarization on the transport properties has been clarified. It is 
shown that the retro-refiectivity of the standard Andreev reflection process is broken in the 
presence of an exchange field and that the surface bound states due to superconducting pair 
potentials do not contribute to the spin current. Next, the ferromagnetic insulator including 
the spin-filtering effect are analyzed. It is shown that the spin-polarization gives asymmetric 
peak splitting. Moreover, various features in the splitting of ZBCP due to the ferromagnetic 
insulator, the Zeeman splitting, and the BTRS states effects are analyzed in detail. It is 
shown that the spin-polarized tunneling gives quite important information to identify the 
origin of the ZBCP splitting. By comparing the present analysis with experimental data, 
we expect that the mechanism of the peak splitting in high-Tc superconductors will be well 
identified. In the present model, we have neglected the effects of spin-orbit scattering^ 
and the non-equilibrium properties of superconductorsH'ii. Inclusion of these effects would 
be necessary for a complete theory. The formulation for triplet superconductors will be 
presented in another publicationS. 
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FIGURES 



FIG. 1. Schematic illustration of the elastic reflection of quasiparticles in the F/I/S junc- 
tion. For all trajectories, momenta parallel to the interface arc conserved. This means that the 
retro-reflection property of Andreev process is lost due to the exchange interaction. In the flgure, 
the anisotropic pair potential of d^2_y2-wave symmetry is also shown. 

FIG. 2. The normalized conductance spectra for the charge current (Tq{E) as the function of 
X = U/EpN with (3 = and Zq = (the transparent limit). As X becomes larger, the peak 
around zero-bias level is largely suppressed. 

FIG. 3. The normalized conductance spectra for the charge current (Tq{E) as the function of 
X with P = tt/4 and Zo = 5. As X becomes larger, the height of the ZBCP is largely reduced. 

FIG. 4. The comparison between the normalized conductance spectrum for the charge current 
<7q{E) and that for the spin current (Js{E) for X = 0.7, (3 = 7r/4 and Zq = 5. Since the ZBCP 
originates from the current carried by the surface bound states, the peak disappears for the spin 
current. 

FIG. 5. The normalized conductance spectra for the spin current (Ts{E) as the function of X 
with /? = and Zq = 5. As X becomes larger, the spin current is increased. Note that the peak at 
£^ = is larger than unity when X is close to one. 

FIG. 6. Orientational dependencies of Jp{E) for d^2_y2-wave superconductors for E = 2Aq 
and X = 0.7 are plotted for various Zq values. Closed dots in the figures correspond to those for 
s-wave superconductors. The large deviations of d^2_y2-wscve from s-wave for small values of Zq 
are originated from the distribution of the pair amplitude in fc-space. 

FIG. 7. The effects of ferromagnetic insulator on the charge current for X = and (3 = tt/A. 
When Zq^i = 5 and Zq^^ = 5, a large ZBCP exists. The difference in Zq^-^ and Zq^i induces the 
peak splitting (.^o.t = 3 and Zq^i = 7). As the difference becomes larger, the ZBCP split into 
two peaks and the amplitude of the splitting becomes larger and the peaks become smaller and 
broader {Zq^^ = 2 and Zq^^ = 8). Finally, the peaks in the gap disappear as the difference becomes 
prominent (.^o,T — ^ Zq.j. = 10)- 

FIG. 8. Four types of normalized conductance for Zq^^ = 2, Zq^^ = 8, X = and (3 = 7r/4. 
The normalized charge current (Jq{E) has splitted peaks. The peak of the lower energy and higher 
energy are originated from the up spin component (jq^^{E) and down spin component (Tg^j^(-E'). 
These peaks do not appear in the spin current conductance as{E). 
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FIG. 9. The normalized conductance spectra CTq{E) as a function of X for Zq.^ = 2.5, Zq^i = 7.5 
and /3 = 7r/4. As X becomes larger, the higher energy peak becomes smaller. Thus, the polarization 
can be estimated from the ratio of two peak heights. 

FIG. 10. The normalized conductance spectra crq(E) in an applied magnetic field 
{^gH/Ao = 0.15) as a function of X for Zq^^ = 5, Zq^^ = 5, and (3 = tt/A. As X becomes 
larger, the peak with higher energy is largely reduced. 

FIG. 11. The normalized conductance spectra crq{E) for the BTRS states (A^/Ao = 0.15) as 
a function of X for Zo.j = 5, Zq^i = 5 and /3 = 7r/4. As X becomes larger, the heights of the two 
peaks are reduced symmetrically. 
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